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Abstract. This paper investigates the physical reasons for the 
existence, the tightness and the near universality of the FIR/ra- 
dio correlation for late-type field galaxies whose emission is not 
dominated by an active nucleus. We develop theoretical models 
for the radio and far-infrared (FIR) emission of such normal 
galaxies and study the influence of their main parameters on 
the ratio of the two emissions. In addition, data are used from a 
sample of 1 14 late-type galaxies which allow an estimate of the 
mean energy densities of the radiation field and the magnetic 
field, the latter crudely calculated from the synchrotron lumi- 
nosity using the minimum energy condition, and of the dust 
opacity. These data reveal, for the first time, a reasonably good, 
linear correlation between the energy density of the radiation 
field and the energy density of the magnetic field. Interestingly 
this implies that the two most important energy loss rates for 
electrons, synchrotron and Inverse Compton losses, are pro- 
portional to each other. As a consequence the radio synchrotron 
emission is proportional to the total flux of radiative energy loss 
from the nonthermal Cosmic Ray electrons of a given energy. 
Furthermore we find that on average the galaxies are marginal- 
ly optically thick for the non-ionizing UV light. Including their 
extended magnetic halos, galaxies are also found to be on av- 
erage marginally optically thick regarding the radiative energy 
losses of the radio synchrotron emitting, non-thermal Cosmic 
Ray electrons. Exceptions may be galaxies with a very low 
(compared to our Galaxy) present star formation rate. Combin- 
ing these semi-empirical results with the theoretical emission 
models we show first of all that in the optically thick case the 
linear correlation between the energy densities of the radiation 
field and the magnetic field is a necessary condition for the 
existence of the observed FIR/radio correlation. Secondly we 
show that the individual dispersion of the FIR to radio con- 
tinuum ratio caused by uncertainties of any single one of the 
parameters considered is significantly less than the empirical 
dispersion of the FIR/radio correlation. The radio and the FIR 
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emission are mainly determined by their sources, i.e. massive 
Supernova precursor stars, because galaxies act in a reasonable 
approximation as calorimeters for the stellar UV radiation and 
for the energy flux of the Cosmic Ray electrons produced in 
their disks. This constitutes a very general theoretical explana- 
tion for the FIR/radio correlation. In principle one can now go 
a significant step further and use this observed correlation plus 
the theory and the parameters of one galaxy as a normaliza- 
tion to deduce in turn the mean magnetic field strength of other 
normal galaxies simply from their total mean radiation energy 
density. 

Key words: galaxies - radio continuum - far-infrared - star 
formation - galaxies: magnetic field 



1. Introduction 

The nearly universal correlation between the integrated far-in- 
frared (FIR) and radio continuum emission of late-type field 
galaxies (see Helou 1991 for a review of the observations) is 
best considered in the framework of a common dependence of 
the FIR and the radio continuum emission on the formation 
of massive stars: (i) the FIR emission is mainly due to dust 
heated by the ionizing and non-ionizing UV radiation from 
these stars, (ii) their ionizing radiation is the source of the 
thermal radio emission, (iii) the Supernova (SN) explosions 
which occur at the end of their lives accelerate the Cosmic Ray 
(CR) electrons responsible for the nonthermal radio emission 
(see Condon 1992, and Volk & Xu 1994, for reviews). However, 
the actual existence of a correlation and especially its extreme 
tightness (~ 0.2 in the logarithm, Wunderlich et al. 1987) which 
is the strongest among various correlations found for many 
different radiations has been difficult to understand, given the 
very different radiation mechanisms of the FIR and the radio 
continuum. 

So far the problem has been addressed, in a qualitative man- 
ner, in two different asymptotic theories assuming exactly the 
opposite extreme conditions. In an 'optically-thick' picture, the 
so-called "calorimeter theory" (Volk 1989; Xu 1990; Lisenfeld 
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tence and the tightness of the correlation may be explained by 
three plausible but far reaching properties of galaxies: First of 
all the main part of the relevant stellar radiation, presumably the 
UV radiation from massive stars, is absorbed by the dust grains 
within a galaxy. Secondly, also the energetic CR electrons that 
are produced by eventual SN explosions of these stars lose their 
energy mainly within the galaxy and its halo, predominantly by 
Inverse Compton collisions with the photons of the interstellar 
radiation field as well as by synchrotron radiation in the magnet- 
ic field. In this case the galaxy as a whole acts like a calorimeter 
for its own primary photon and electron emission, independent 
of the details of production and transport of these two carriers 
of energy. Finally, the energy density of the interstellar radia- 
tion field, U ra d, is in a constant ratio with the magnetic field 
energy density, U B . This means that the radio synchrotron e- 
mission is proportional to the total (radiative) energy loss rate, 
and that both FIR and nonthermal radio emission are essential- 
ly proportional to the formation rate of the massive stars with 
but a weak dependence on all other parameters. The result is 
a tight, roughly linear correlation. When the contribution from 
the old stars to the heating of dust is taken into account, the 
weak nonlinearity of the correlation (Pi. 49GHz oc L^,) can also 
be explained satisfactorily (Xu et al. 1994). 

Alternatively, in the 'optically-thin' picture (Helou & Bicay 
1993; Chi & Wolfendale 1990), it is assumed that most of 
the radiation which is responsible for the heating of dust, and 
most of the CR electrons which produce synchrotron radiation, 
escape from a typical galaxy disk, i.e. disk galaxies are optically 
thin for both the primary photons and CR electrons. If this 
is indeed the case for the majority of late-type galaxies, it is 
argued (Helou & Bicay 1993) that the FIR/radio correlation 
indicates a proportionality between the two escape rates, namely 
to oc t esc /t sync where to is the optical depth of the disk for 
the primary photons, and t esc and t sync the escape time scale and 
synchrotron loss time scale of CR electrons. The nonlinearity of 
the FIR/radio correlation is then interpreted as a consequence 
of the dependence of the CR electron escape rate on the FIR 
luminosity (Chi & Wolfendale 1990). 

In this paper we shall study the FIR/radio correlation with 
a quantitative perspective and scrutinize these asymptotic the- 
ories. Theoretical models for both the radio continuum and the 
FIR emission will be developed. We will show explicitly that 
a large number of parameters are involved. This demonstrates 
the significance of the correlation. For the major parameters 
determining the emissions, e.g. the optical depth, the magnetic 
field strength, the ratio between the energy densities of the ra- 
diation field and the magnetic field, etc., we will avoid a priori 
assumptions but rather constrain their mean values and disper- 
sions using observational data from a sample of 114 nearby 
spiral galaxies and other data from the literature. It is clear that 
such a correlation can only exist for normal field galaxies, out- 
side the core of galaxy clusters, because the action of an active 
nucleus (AGNs) or the intense interaction with other galaxies 
would introduce effects independent of the star formation and 
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an AGN as well as clustered galaxies. 

In Sect. 2 we present the models for the FIR and the radio 
emission. In Sect. 3 we discuss observational constraints for 
some of the parameters, especially the ratio U ra d/UB ■ In Sect. 4 
the sources of the dispersion of the FIR/radio ratio are evaluated, 
whereas Sect. 5 presents our results concerning the question of 
the optical depth of galaxies. Sect. 6 contains the discussion. In 
Sect. 7 the main conclusions are summarized. 

2. Models of the radio continuum and the FIR emissions 

Many parameters, most of which do not depend directly on 
the formation rate of massive stars, are involved in the radia- 
tion processes of the radio continuum and of the FIR emission 
of galaxies. In order to study their influence on the FIR/radio 
flux ratio, we give in the following a detailed theoretical de- 
scription of the mechanisms of the two emissions, and evaluate 
the importance of individual parameters for the corresponding 
process. 

2.1. Dust heating and FIR emission 

The FIR radiation of normal galaxies is completely dominated 
by the thermal emission of dust heated by stellar radiation. The 
FIR luminosity (in the wavelength range of 40-120/tm) of such 
a galaxy can be derived from the following formula: 

/>m2 />t'(m) 

Lfjs. = / C(m)$(m) dm / r(t) dt 

Jml JO (J) 

J P\(m) f x (m) d\ 

with m\ = 1M© and TO2 = 100M Q specifying the mass range 
of stars which are responsible for the heating of dust (we ignore 
stars less massive than the Sun because they contribute little to 
the dust heating, and take an upper cut-off of 100M Q for the 
Initial Mass Function (IMF) $(m)); t'(m) is the main sequence 
life time of a star of mass m, /a (to) the radiation spectrum of 
a star, and ((m) is the fraction of the absorbed radiation from 
a certain star to be re-radiated within the wavelength range 
40-120/xm. From equation (1), the FIR luminosity depends in 
addition on the star formation rate and its time dependence r(t), 
where t is the looking back time, as well as on the probability, 
P\(m), of radiation at wavelength A, radiated by the star of 
mass m, to be absorbed by dust. 

In order to parameterize these dependences, we make the fol- 
lowing assumptions: 

(1) The IMF is a power law between mi = 1 M Q andTO2 =100 
M Q , with power index k : 

$(m) = $ m~ k (1 < m/M Q < 100). (2) 

(2) The stars in the mass range 1 — 100 M Q can be di- 
vided into 3 populations: i) the massive ionizing stars of 
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so that Jq <m) r(t) dt ~ ro x i'(m), where ro = r(i = 0) is 
the current star formation rate, ii) the intermediate massive 
stars of 5 < to < 20M Q which are mostly responsible for 
the non-ionizing UV radiation of a galaxy, with a charac- 
teristic life time t'(m) = 10 8 years, hi) the low mass stars 
of 1 < m < 5M which are mainly responsible for the 
blue radiation of a galaxy, with a characteristic life time of 
t'(ni) = 3 10 9 years. For populations ii) and iii) the time 
integration of the star-formation rate can be approximated 
by r$ x t'(m) and rg x t'(m), where rs and are the star 
formation rate averaged over the last 10 s years and 3 10 9 
years, respectively (Xu et al. 1994). 
(3) The massive ionizing stars heat the dust locally, i.e. they 
are the sources of the HH-region-associated FIR emission. 
The absorption probability of their radiation is taken as a 
constant (Xu et al. 1994): 



P x (m) = P HII = 0.6 



(m > 20M Q ) 



(3) 



Using the model of Desert et al. (1990) for the dust heated 
by an 05 star (dilution factor X=0.01), we estimate that the 
fraction of the HH-region-associated dust emission through 
the window 40-120/um amounts to 



am) = Chii = 0.5 



(m > 2OM ). 



(4) 



(4) Dust heating due to intermediate massive stars and old s- 
tars occurs through the diffuse and approximately uniform 
interstellar radiation field which heats the diffuse dust in 
the disk. The resulting FIR emission is the so-called Cirrus. 
We take the corresponding ( value in all galaxies as that of 
the cirrus emission in the solar neighbourhood (Desert et al. 
1990), which is 30 %: 



C(m) = Cdiff = 0.3 



(to < 20M Q ). 



(5) 



(5) The radiative transfer for the non-ionizing emission of the 
intermediate massive stars and the old stars is modeled using 
an infinite parallel-slab geometry, assuming that stars and 
dust are uniformly mixed. According to Xu and De Zotti 
(1989) the absorption probability is then estimated from the 
optical depth of galactic disks: 



P\(m) = 



0.6 + 



r /0.9 



(m < 2OM ) 



(6) 



where t' x = (3t\ is the so called 'effective optical depth' of 
the disk which takes into account the effect of scattering, 
and 1 — a\ < (3 < 1, where a\ ~ 0.6 is the albedo (Mathis 
et al. 1983). We take /3 = 0.7. For the sake of simplicity, we 
make the further approximations: 



P\(m) ■■ 



/0.9 

u.u-r / 
/0.9 



B 



0.6+r 



7TCT 



(5 < to < 20M Q ) 
(1 < to < 5Mq) 



(7) 



where t' uv is the effective UV optical depth at 2000A, 
and t' b is the effective blue optical depth at 4300A. The 
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adopted extinction curve (Savage and Mathis 1979). This 



gives t' b = 0.47 r ( 



uv 



With the above assumptions and simplifications, equation (1) 
can be rewritten as: 



£fir =ro Phh Chii Liik) + r 8 Puv(ruv) Cdiff Lu{k) 
+ r 9 P b (t b ) Cdiff Liu{k) 

where 



Li(k) = / <&Qm~ k t'(m)L(m) dm 

J2QM 



(8) 



(9) 



is the total luminosity of massive ionizing stars for ro = 1, 
and L(m) = J f\(m) dX is the total luminosity of a star of 
mass m. Since the product of the life time and the luminosity 
for stars of a given mass is well determined from the theory of 
stellar evolution, L, is only a function of the parameter k which 
specifies the IMF. The same is true for the similar variables Lu 
(for the intermediate massive stars - between 5 and 20 M Q ) 
and Lm (for the old stars - between 1 and 5 M ). According to 
Xu et al. (1994) we will assume for the total radiation emitted 
during a stellar main-sequence life-time: 



t'L = 



_ / 10 9 95 to L & yr 



10' 



9.6„ 



» 3 /2 L 



(1M < m < 5M ) 
(5M q < m < 100M G ) 



(10) 



In Lisenfeld (1993) a more detailed model for the dust 
heating has been considered in which the dust absorption of 
non-ionizing radiation (A > 912A) and the dust absorption of 
ionizing radiation in HII regions are explicitly taken into ac- 
count. The fraction of stellar radiation absorbed by dust could 
therefore be calculated without the need of simplifying assump- 
tions like Eqs. (3) and (7). However, the results of Lisenfeld 
(1993) support the validity of these assumptions: Taking as an 
example tjjv = 1> Eqs. (3) and (7)) predict that in the mass 
range of 1-5 M Q (respectively 5-20 M Q , 20-100 M Q ) a fraction 
of 37 % (54 %, 60%) of the radiation of a star is absorbed. 
The more detailed calculations yield for the same mass ranges 
25 -40 % (40-60 %, 60-65%), where the smaller numbers of 
theses ranges refer to the low mass ends of the mass intervals, 
and the larger numbers to the high mass ends. For different dust 
opacities the agreement is equally good. 

2.2. Radio continuum emission 

2.2.1. Synchrotron emission of CR electrons 

In a one-dimensional, steady state diffusion model, the diffusion 
equation describing the production, diffusion and energy losses 
of electrons reads as follows (e.g. Berezinsky et al. 1990): 



D(E) 



d 2 f(E,z) _ d 



dz 2 



+ — \b(E)f(E,z))=-Q(E,z) , (11) 



where z is the height above the galactic midplane, and f(E, z) 
denotes the phase space density of the electrons, i.e. f(E, z) 
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(E,z). 



D(E) = D 



E 



lGeV 



(12) 



is the energy dependent diffusion coefficient, and b(E, z) de- 
scribes the electron energy losses in the energy range considered 
here (E > 1 GeV). These are dominantly Inverse Compton and 
synchrotron losses: 



dE 



syn+IC 



= -a T c 



E 

m e c 2 



(U md + U B ). (13) 



Here, er T is the Thompson scattering cross section, c is the speed 
of light; Ub denotes the energy density of the magnetic field and 
[7 ra d the energy density of the radiation field below the Klein- 
Nishina limit (e.g. Longair 1992). For the electron energies 
considered here (E > 1 GeV), this limit for the photon energies 
lies in the X-ray domain. Thus, to a very good approximation 
[/rad represents the energy density of the entire radiation field. In 
order to simplify Eq. (1 1), we make the following assumptions: 
1) The CR electrons are produced only by Supernova Rem- 
nants (SNRs) from massive stars which are situated in the 
galactic midplane at z=0, and the production rate is propor- 
tional to the SN rate. The source function has a power-law 
spectrum (e.g. Volk et al. 1988): 



Q(E,z) = 8(z) 



E 



m e c 



(14) 



where z^n is the SN rate, </sn is the average of the number of 
relativistic electrons that are produced per energy interval 
dE by a SN, and x > 2 denotes the spectral index. We have 
therefore neglected other possible CR electron production 
mechanisms, e.g. shocks associated with type I SN from 
(old) White Dwarfs in accreting binary systems (Xu. et al. 
1994), stellar winds (Cesarsky and Montmerle 1983), spiral 
arms (Duric 1986), and pulsars (e.g. Volk 1992). Assuming 
that stars more massive than 5M are precursors of the SN 
considered, v$n is determined by the adopted IMF (i.e. 
the parameter k) and by r%, the star formation rate averaged 
over the last 10 8 years, (the life time of intermediate massive 
stars which considerably outnumber the massive stars): 
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nNsN(k) 



1) is the number of SN 



where NsN(k) 
when r(t) = 1 yr~ 
2) The energy densities C/ ra d of the radiation field and Ub of 
the magnetic field (and therefore also the field strength B) 
are step functions of the height z over the galactic plane: 
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Fig. 1. Geometry of the magnetic field B and the energy density of the 
radiation field [Trad in the galactic halo: The magnetic field B a and the 
energy density of the radiation field (?7 ra d)<i in the outer halo (z > L) 
are by a factor S and p, respectively, lower than the values in the inner 
halo, Bi and (frad)t- 
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(\z\ > L), 



(16) 



where L (~ 1 kpc) is a free parameter to separate the so 
defined outer halo from the inner halo, and 6 < 1 as well as 
p < 1 (Fig. 1). This is a first order approximation to describe 
the actual spatial decrease of the two fields with distance 
from the galactic disk. 
3) A free escape boundary exists at | z \ = z h , the outer boundary 
of the halo: 



f(E,\z h \) = 



(17) 



(e.g. Ginzburg & Syrovatskii, 1964, p. 300). An approxi- 
mate solution of Eq. (1 1) for a uniform halo (i.e. S = p = 1) is 
given by: 
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(18) 



with 
P=~ 



2x + p — 3 
2(1 - p) 



(19) 



iF\(a,b,x) denoting the confluent hypergeometric function (or 
Kummer function) (e.g. Mathews and Walker, 1970), and 
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These variables are proportional to the ratio of the diffusion 
time scale Tnjf to the energy loss time scale ti oss : 
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and 
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(25) 



For a nonuniform halo, a similar approximate solution can be 
found with an analogous substitution by solving Eq.(ll) sep- 
arately for the inner and outer part of the halo. The solutions 
are matched at z = L such that both the density f(E, z) and the 
flux D(E)(df(E, z)/dz) are continuous (Lisenfeld 1993). The 
accuracy of these solutions is discussed in Appendix A. 

The total synchrotron emission P syn (v) of a galaxy is cal- 
culated as: 



-I 



Zh dE 
dzf(E,z) ' 



dt 



dE 



(26) 



Here we use the simplifying assumption that an electron emits 
its synchrotron radiation completely at the maximum frequency 
v of its synchrotron spectrum: 



v 



1 



where v G = eB/m e c 2 is the gyrofrequency. 

We assume a constant spectral index of the injection spec- 
trum x = 2.2 (Volk et al. 1988) although theoretically x = 2.1 
(Berezkho et al. 1993) is an equally plausible value (which then 
requires \x ~ 0.6 to establish agreement with CR observations 
in our Galaxy). 

We further assume that each SN (mainly type II and lb) 
is a "standard candle", i.e. they all produce the same amount 
and the same spectrum of relativistic electrons (Markiewicz, 
Drury & Volk 1990). Although this might not hold precisely for 
individual SN, it is probably a good assumption when averaged 
over a whole galaxy since at any given epoch ~ 1000 SNRs act 
as CR sources taking our Galaxy as a typical case (with an active 
life time ~ 10 4 to 10 5 yrs for a SNR and v m ~ 1 /30yr). Under 
this assumption, the parameter qgN is a constant. Its value 
can be deduced from the data of the Galaxy: The intensity of 
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W/Hz (extrapolating the results of Beuermann et al. 1985 and 
assuming a nonthermal spectral index of 0.9), the SN rate is, 
within a factor of about 2, vgN = 0.04 yr _1 (van den Bergh 
and Tammann, 1991), the strength of magnetic field is B = 
5/jG (Sofue, Fujimoto & Wielebinski 1986) and Do = 10 29 
cm 2 /sec, Zh = lOkpc (Berezinsky et al. 1990). From these data 
we derive the value of q SN = 1.6 10 48 eV -1 . The value of q S N 
is, mainly because of the uncertainty in the Galactic SN-rate, at 
least uncertain to within a factor of 2 (0.3 in the logarithm). 

We can therefore parameterize the total synchrotron emis- 
sion of a galaxy in the following symbolic form: 

Psyn(v) = QSN N 

TlAx, {Ur^/U^iiBi, D , z h , S, p, L), 

where n„ is a nonlinear function of the nine parameters indi- 
cated. 

It is instructive to consider P syn in the asymptotic limits 
\sh\ ^> land|s/j| <C 1 . This requires the asymptotic expansions 
of the confluent hypergeometric function: 
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where T{x) denotes the gamma function (e.g. Mathews & Wal- 
ker, 1970). 

For a homogeneous halo we obtain for \sh\ <C 1: 



Psyn(v) = —77; Cl (m e C ) Z h \ 

2D v G \v G J 



(29) 
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\sh(E)\ -C 1, means that the electrons have not suffered any 
substantial energy losses when they reach z h . This is reflected 
in Eq.(29) by the fact that the synchrotron emission does not 
depend on the energy losses. The diffusion coefficient and the 
halo size z h , however, are important parameters in this 'optically 
thin' limiting case. 

In the opposite case, |s ft | > 1, the synchrotron emission is 
given by: 



Psyn = I'SNQSN ' 



(m e c z ) 



.2\2 



4(x -\)v G V v G ) 1 + U rdd /U B 
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oc 



(30) 



In this case of high frequency v, energy losses dominate the 
electron propagation. The electrons loose their energy down to 
E or below before they reach the outer boundary of the halo. 
The strength and energy dependence of the diffusion coeffi- 
cient are not important for the total synchrotron emission (for 
a homogeneous halo) and therefore neither Do nor \i appears. 
This asymptotic case corresponds to the "calorimeter" theory 
considered by Volk (1989). It is characterized by a weak de- 
pendence on the magnetic field. In addition, P syn is inversely 
proportional to (1 + [/ ra d/%)- 
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Especially at high frequencies, thermal bremsstrahlung of ion- 
ized gas contributes to the radio emission as well. It is mainly 
determined by the number of Lyman continuum photons, Niy C 
(in sec -1 ), emitted by exciting stars: 
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(31) 



[W/Hz] 



(e.g. Condon, 1992), where f ion is the fraction of Lyman contin- 
uum photons absorbed by the gas, T e is the electron temperature 
and TVLyc is a function of the IMF as well as of the present star 
formation rate ro: 
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20M Q 

ron Lyc (k) 



$om nL yc (m)s(t = t'(m)) dm 



(32) 



where fiL yc (m) is the production rate of Lyman continuum 
photons by a star of mass m, and n,L yc = $o f nL yc (m)m~ k 
is the total production rate of Lyman continuum photons when 
r(t) = 1 yr~ l . The values of nL yc (m) are taken from Giisten 
& Mezger (1983). For the other parameters of the thermal 
bremsstrahlung we will assume standard numbers: 

1) T e = 1000K, a typical temperature for the electron gas in a 
HII region (e.g. Osterbrok, 1989) 

2) fion = 0.8 (Mezger 1978). 

We can therefore parameterize the total thermal radio emis- 
sion of a galaxy as follows: 

PthM = ron Lyc (k)q th (f ion , T e ,v) , (33) 
where q th is a known nonlinear function of its arguments. 

2.3. The FIR/radio ratio 

From sections 2.1 and 2.2, the FIR luminosity Lfir and the 
radio continuum power P v , and consequently the logarithm of 
their ratio q v = \og(L F i R /P v ), can finally be parameterized as 
follows: 

LmL =ro Phii (hii L^k) + r 8 Puv(tuv) (diff Lu(k) 

+ r 9 Pb(t b ) (diff Lui(k) (34) 

and 

P v =Pth(v) + P syn (v) 

=roq t h(fion, T e ,v)n Lyc (k) + r s N S N(k) qsN (35) 
U v (x, p, (U rd d/U B )i, Bi, D , z h , S, p, L) 

In total there are 20 parameters involved in q v : 
1) Parameters concerning both emissions: 
- the star formation history: ro jr% 



2) Parameters concerning the FIR emission: 

- dust absorption: tb , Wv , Phii 

- dust emission: (diff, (hii 

- star formation history: rg j r% 

3) Parameters concerning the radio emission: 



thermal bremsstrahlung: 
electron production: 
energy losses: 
electron diffusion: 
geometry of halo: 



fion ; T-e 

x,qsN 

(U rd d/UB)i, Bi 

Do,n,z h 

L,p,6 



For some of these parameters, assumptions concerning their 
values have been made: In section 2.1 we have assumed that 
P HII = 0.6, ( HII = 0.5 and ( diff = 0.3, t b = 0.47 x t uv , and 
in section 2.2 we have assumed that x = 2.2 (or 2.1), qsN = 
1.6 10 48 eV" 1 , T e = 7000 K, and f ion = 0.8. These parameters 
should be more or less the same for galaxies, and they affect 
the FIR/radio ratio only weakly. For example, the parameters 
Phii, T e , and f ion depend only on the local conditions for 
HII regions, and affect only the thermal radio emission and the 
'warm' FIR component which are minor components in the 
radio continuum and FIR emissions. We will therefore in the 
following take their values as fixed and neglect the influence of 
their variations on the FIR/radio correlation. 

The absolute values of the star formation rate, i.e. ro,rg 
and rg, appear linearly in both expressions of Lfir and P v . 
Thus q v depends only on the star formation history, i.e. on the 
ratios ro/rg and rg/rg. We will furthermore assume an average 
value ro/r% = 1, corresponding to galaxies with a constant 
SFR during the past 10 8 years. We therefore exclude from this 
analysis starburst galaxies that have had a strongly increased 
SFR during the last ss 10 7 years. 

With these assumptions, there are 1 1 parameters left for a 
more detailed discussion. They are: k, tuv, rg/n, {U rd &lU%)i, 
Bi, Dq, [i, Zh, L, p, and 6. If either one of the asymptotic cases, 
i.e. the 'optically-thick' or 'optically-thin' case discussed in 
the introduction (Sect. 1), applies for most of the late-type 
galaxies which follow the FIR/radio correlation, the number of 
the parameters involved can be further reduced: 
(1) Optically-thick case: In this asymptotic case tuv » 1 
and \sh\ >> 1, i.e. galaxies are optically thick for stellar 
radiation and synchrotron electrons lose most of their energy 
before escaping from a galactic halo. Consequently, P uv ~ 
Pb — 1 in Eq.(34), and the second term on the R.H.S. 
of Eq.(35), i.e. the synchrotron radiation power P syn , is 
approximated by expression (30). The FIR/radio ratio is 
therefore a function of only 4 parameters, namely the slope 
of the IMF k, the star formation history parameter r% jrg, the 
strength of the magnetic field B, and the ratio {Uad/U^i- If 
we further ignore the contribution of old stars to the heating 
of dust (the third term on the R.H.S of Eq.(34)), then the 
model reduces to the 'calorimeter theory' given by (Volk 
1989), where the FIR/radio ratio is inversely proportional 
to (1 + U la d/UB)i, and depends only very weakly on the 
magnetic field strength B and the parameter k. 
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Then, py v and oc T^y- (c.f. Eq.(7)), and P s?/ „ is approx- 
imated by Eq.(29). The FIR/radio ratio is then a function of 
9 parameters: the UV optical depth iyv> me magnetic field 
strength B and parameter S, the four parameters describing 
the electron propagation: Dq, z h , p, L, and the parameters 
k and rg/rg. 



3. Observational constraints on the parameters 

3.1. Parameters concerning both emissions 

1) For the relevant mass range 1 < m/M Q < 100 we adopt 
a slope of the IMF of k = 2.5, consistent with the IMF 
determined for the solar neighbourhood by Miller & Sca- 
lo (1979). Scalo (1986) compared the IMF's of different 
galaxies and concluded that there is no observational evi- 
dence for the slope of the IMF to change between galaxies 
by more than ±0.5. 

3.2. Parameters concerning the FIR emission 

2) Assuming a constant gas-to-dust ratio with the same value 
as observed in our Galaxy, we are able to give an estimate 
for the average dust opacity < t uv > and its dispersion for 
our galaxy sample. Details are described in Appendix B. 
The result is < tuv >= 1 • 1 with a logarithmic dispersion 
of6i ogr = 0.27. 

3) Following the argument of Xu et al. (1994) we assume that 
the average value of rg/r& = 1. Using the formulae given 
in Xu et al. (1994), we find a dispersion of the logarithm of 
this ratio of 0.45. 

3.3. Parameters concerning the radio emission 

4) The average values and dispersions of the magnetic field and 
of the energy density of the radiation field can be estimated 
for the galaxy sample. The radiation field of a galaxy can be 
calculated from its bolometric luminosity L\> i and its radius 
R. For the magnetic field, only a rough estimate is possible 
from the total nonthermal radio emission through the min- 
imum energy assumption. The detailed procedure of the 
calculation of the two quantities is described in Appendix 
C. For our galaxy sample, the average resulting energy den- 
sity of the radiation field is (U ra d) = 1.4eV/cm 3 and the 
average value of the magnetic field is (B min ) = 1.5pG with 
a dispersion in the logarithm of Si og B = 0.12. 

In Fig. 2 the energy densities of the magnetic field and the 
radiation field are plotted. Clearly a correlation can be seen. 
A linear regression analysis yields: 



C/ B ex U^l 



28±0.12 



(36) 



with a correlation coefficient of r = 0.74. The average ratio 
of the two quantities is 
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Fig. 2. Energy densities of the radiation field and energy densities of 
the magnetic field, for the galaxy sample 



The dispersion in the logarithm of this ratio is S\ og x =0.16. 
Thus we conclude that the two quantities are reasonably 
well correlated with a roughly linear relation and that their 
values are nearly equal. The slight nonlinearity in the re- 
lation will not be studied further because it would be an 
overinterpretation of the simple approach used to calculate 
U rd d (see Appendix C for details). 

5) The parameters concerning the electron diffusion (Do, p, 
Zh) have been estimated from observations for our Galaxy. 
We will adopt these numbers: D(lGeV) = 10 29 cm 2 /s and 
Zh = 10 kpc (Berezinsky et al. 1990), and the energy depen- 
dence of the diffusion coefficient p = 0.5 (e.g. Berezinsky 
et al. 1990). This is consistent with a value x = 2.2 of the 
source spectral index (see section 2.2). 

6) The extension of the inner halo above the galactic disk is 
assumed to be L = \kpc. This value lies within the range 
determined for the extension of the thin radio disk for the 
edge-on galaxies NGC 819 (L=400 pc, Dahlem, Dettmar & 
Hummel 1994) and NGC 4631 (L=1.5 kpc, estimated from 
the maps by Golla & Hummel 1994). 

The parameter p, describing the decrease of the radiation 
energy density away from the galactic disk is given by the 
geometrical dilution. With the above values of L and z h , we 
derive p = 0.4. 

The decrease of the magnetic field in the halo is more diffi- 
cult to estimate. For our Galaxy, dynamo theory predicts a 
ratio of the magnetic field strength in the halo, Bh, to that 
in the disk, B d , of B h /B d « 0.1 (Ruzmaikin et al. 1988) 
In other galaxies this ratio might be higher: For NGC 463 1 
a ratio of Bh/Bd « 1 is conceivable (Hummel, Beck & 
Dahlem 1991). We will consider various values for S in this 
range. 
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Fig. 3. Spectral index a ayn of the synchrotron radiation shown for S = 
1 , p = 1 (full line), 5 = 0.5, p = 0.4 (dashed line) and 5 = 0.2, p = 0.4 
(dashed-dotted line). B{ = IpG, {U r a&/Un)i = 1.1, according to the 
results of Sec. 2.2 and Zh = lOkpc, Do = 10 29 cm 2 /s, according to the 
values in our Galaxy, and L = lkpc are adopted. 

A further constraint to the parameters concerning the syn- 
chrotron emission, is given by the nonthermal spectral index, 
a syn . Klein (1988) separated for 13 galaxies the thermal from 
the nonthermal radio emission and derived an average non- 
thermal spectral index of 0.88 ± 0.06. Independently, Duric, 
Bourneuf & Gregory (1988) obtained an essentially identical 
result. 

The prediction for the spectral index within this model is 
given in Fig. 3 for the parameters listed above. In addition to 
S = 0.2, p = 0.4 and S = 0.5, p = 0.4 we have plotted for com- 
parison the case of a homogeneous halo and disk (6 = p = 1). 
The predicted spectral index of the synchrotron radiation be- 
tween 100 MHz and 10 GHz lies between 1.0 — 1.1 in the uni- 
form case, and 0.9 - 1 .0 for an outwardly decreasing magnetic 
field. With the slightly lower source spectral index (x = 2.1) 
advocated by Berezkho et al. (1993), these values decrease to 
0.95 - 1.05 in the homogeneous case, and to 0.85 - 0.95 for a 
decreasing magnetic field. The decrease of the magnetic field in 
the halo produces a lower spectral index which is in reasonable 
agreement with the observations, indicating that the parameter 
values we have chosen are in the right range. 

4. Dispersion of the FIR/radio correlation within the model 

For our galaxy sample we find a mean of the logarithm of 
the FIR/radio ratio q\A9GHz = 15.06 with a dispersion of 
a qiA9GHz = 0.24. This can be compared to the dispersion de- 
rived within the model: We calculated the derivative of q\A9GHz 
with respect to the logarithms of various parameters, represent- 
ed by Xi (i=l, 2, ...), using the mean values of the parameters 
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quantitative estimate how sensitive the FIR/radio ratio depends 
on the respective parameter. For some of the parameters we 
could derive the mean value and its dispersion, 6\ og ( Xi ), and thus 
we can give an estimate for their contribution, a Xi , to the scatter 
in the logarithm of the FIR/radio ratio Q\.a9GHz by: 

d \og{xd 

where the derivative is taken at the radio frequency of 1 .49 GHz. 
Table 1 summarizes the results (for S = 0.2, p = 0.4). The pa- 
rameter k is not included because its contribution to the disper- 
sion of the correlation is not significant (<7& < 0.1 for 8u < 0.5). 
The parameters S, p, L and p, are also dropped because their in- 
fluence on the synchrotron emission is only second-order and 
therefore less important compared to (Uai/U^i, Bi, Zh and 
D . 

Table 1. Contributions to the dispersion of the FIR/radio correlation 
by individual parameters 





< Xi > 


a«1.49Gffi 

a loud j ) 






TUV 


1.1 


0.4 


0.27 


0.11 


Bi 


7.5/iG 


-0.9 


0.12 


-0.11 


{Utni/U^i 


1.1 


0.3 


0.16 


0.05 


rg/rg 


1 


0.3 


0.45 


0.14 


D 


10 29 cm 2 /s 


0.5 






Zh 


lOkpc 


-0.4 







For the first four parameters in Table 1 , namely the UV op- 
tical depth Tuv, the strength of the magnetic field in the inner 
thin disk Bi, the ratio between the energy densities of the ra- 
diation field and that of the magnetic field (U m a/UB)i, and the 
star formation history parameter r^jr%, quantitative estimates 
of their contributions to the scatter of the FIR/radio ratio are p- 
resented. The absolute values are all smaller than the dispersion 
of q\A9GHz - The contribution from the parameter (U m a/UB)i is 
negligibly small, due to the reasonably tight correlation between 
U rd d and Ub (see section 3.3) and a relatively weak dependence 
°f qi.49GHz on this parameter (the logarithmic derivative equals 
to 0.3). The contributions of D and z% are not given, because 
we do not have the data to evaluate their dispersions although, 
according to our calculation, the FIR/radio ratio depends on 
these two parameters rather sensitively (the logarithmic deriva- 
tives are 0.5 and 0.4, respectively). 

In principle the variance, i.e. the square of the dispersion of 
qi. 49GHz is given by: 

^■.«», = a *i + 2 a *i ri 'j ■ (39) 

The second term on the R.H.S. takes into account correlations 
between the parameters Xi and Xj, where r^j is the linear cor- 
relation coefficient between x i and x j . In Table 2 we list the 
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eters tuv, B, and rg/r%, for which we find a significant a in 
Table 1. Significant anticorrelations between tjjv an d rg/rg, 
and between B and rg /r% are found, which result in significant 
contributions (with opposite signs) to the variance of Qia9GHz 
(Table 2). The anticorrelation between the optical depth and 
rgjr% may be understood in terms of the well established corre- 
lation between the gas column density and the present day star 
formation rate (e.g. Buat 1992), given that our tjjv is estimated 
from the HI column density. The physical interpretation of the 
second correlation, i.e. the B vs. rg/r% anticorrelation, may be 
drawn from the dynamo theory for the galactic magnetic field 
(Ko & Parker 1990), which predicts an enhancement in B for 
a more active present day star formation activity (i.e. a lower 
rg/r&). No significant correlation is found between tjjv an d B 
although both are correlated with rg /r%. 

Table 2. Contributions to the variance of the FIR/radio correlation by 
three cross correlations 







2&Xi & Xj 


log-rt/v v.s. logB 


-0.015 


0.0004 


\ogruv v.s. log(r 9 /r 8 ) 


-0.50 


-0.015 


log B v.s. log^/rg) 


-0.62 


0.021 



From the values in Table 1 and Table 2 we can predict 
a lower limit of a* ,„„ > 0.052. This is consistent with 

11.49GHz 

the observed value a 2 „, ,„„„ = (0.24) 2 = 0.058, leaving some 
room for the unknown contribution from the parameters Do and 
Zh- The absolute value of q\A9GHz predicted by our model is 
QiA9GHz = 14.6 ± 0.3, slightly lower than the observed value 
q\A9GHz = 15.06 ± 0.24. Several uncertainties in our model 
may be responsible for this discrepancy: (1) A lower-mass cut- 
off of 8M , instead of the 5M used in our calculation for the 
SN rate vsn, will raise the prediction for q\ A9GHz to 14.9±0.3; 
(2) the error in the derivation of qsN- It should be pointed out, 
however, that our results on the dispersion of Qi.49GHz are not 
at all affected by this discrepancy. 

5. Optically-thick or optically-thin? 

With the parameters as they are determined using our galaxy 
sample, we now can discuss the question whether most late-type 
galaxies, as represented by our sample, are 'optically-thick' or 
'optically-thin' for their primary photons and CR electrons. 

The mean optical UV depth tuv = 1.1 with the logarithmic 
dispersion of 0.27 (a factor of 1 .86 in linear scale) indicates that 
most of galaxies are optically thick, although not extremely so, 
for the UV radiation which is most responsible for the heating of 
dust(Xu 1990;Xu&Buat 1995). Our result is based on the dust- 
to-gas ratio adopted from the value in the solar neighborhood. 
The validity of this assumption has recently been confirmed 
by Xu & Buat (1995) and Buat & Xu (1995) who studied the 
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Fig. 4. The fraction of synchrotron emission lost due to electron escape, 
f l S yn , as a function of frequency, for the same parameters as in Fig. 3 

extinction in late-type galaxies using a radiative transfer model 
which allows for the effect of dust scattering and which takes 
input from the UV, optical and FIR observations of a sample 
of nearby late-type galaxies. Xu & Buat (1995) found a mean 
optical depth for the blue radiation (4400A) t b = 0.6 ± 0.04, 

indicating a r s = 1 .08 ± 0.07, consistent with our result. 

& 2000A 

Within our theory, we can also quantify the importance of 
escape of CR electrons: In Fig. 4 the fraction of synchrotron 
emission lost due to electron escape f s ° s ^ is plotted. 

It is defined as: 



- 1 syn 

with P™® x being the maximum synchrotron emission for given 
parameters L, S and p, assuming an infinite extent of the halo. 
At 1 GHz the fraction of synchrotron emission lost is between 
10 — 25%. The dominant part of the energetically possible 
synchrotron emission is thus observed. This result is confirmed 
by the value for \sh \ that we derive in the model: At 1 GHz we 
obtain (in the spatially homogenous case) a value of \sh \ = 0.69 
which means that the diffusion time scale r di ff is 5.5 times 
longer than the energy loss time scale r; oss (see Eq.(23)). 

Therefore we conclude that the 'optically-thick' scenario 
for the FIR/radio correlation describes the situation for late- 
type galaxies better than the 'optically-thin' scenario. On the 
other hand the extreme conditions assumed in the 'calorimeter 
theory', namely t uv >> 1 and f l s ° y s ^ « 1, are not strictly 
fulfilled, as can also be seen from the logarithmic derivatives 
given in Table 1. In the extreme case, one should have near- 
ly zero logarithmic derivatives with respect to tuv, A)> and 
Zh, and a very small derivative with respect to B (Volk 1989). 
Yet we find that the absolute values of these derivatives are 
all significantly larger than zero (Table 1). Nevertheless, these 
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the 'optically-thin' case, namely 0.9 for tuv (Eq(7)), 1 for D 
and 2 for Zh (Eq(29)), and 1 .85 for B. In the optically thin case, 
the stronger dependence of q v on the various parameters would 
result in a higher dispersion: Taking into account only the dis- 
persion for B,tuv an d rg/rg (thus neglecting Do and z^) we 
would already predict of 49GHz = 0.13, exceeding by far the 
observed dispersion of the FIR/radio correlation. Therefore, as 
shown in the previous section, the tight FIR/radio correlation 
can actually be explained by these weak dependencies, as has 
been argued in the 'optically-thick' scenario, without invoking 
any unknown relations among the parameters unless the con- 
tributions of Do and Zh to the dispersion are very significant. 
Here, the tight correlation between D ra d and U B is very im- 
portant for the low predicted value of of ^qhz > reducing the 
contribution of (J7 ra d/%)i to a very small value. Therefore, the 
proportionality between U m( \ and Ub can be considered as a 
necessary condition for the existence of the tight FIR/radio cor- 
relation. Furthermore, we do not find any significant correlation 
between tuv and B, indicating the absence of a correlation be- 
tween the trapping rate of the primary photons and that of CR 
electrons as required by the 'optically-thin' scenario. 

6. Discussion 

6.1. Model simplifications 

In the model describing the radio and the FIR emission, several 
simplifications were introduced, the major ones being the use 
of a 1 -dimensional diffusion model to describe the electron 
propagation and the assumption of an infinite slab geometry for 
the galactic disk in which stars and dust are uniformly mixed. 
In the following we will discuss these assumptions. 

By using a slab geometry for the model of the dust heat- 
ing we assume that the scale heights for dust and stars are the 
same. This is most likely a good approximation for massive 
stars responsible for the UV-radiation (m > 5M©), whereas 
the scale height of lower mass stars is larger than that of the 
dust. Therefore the FIR rises less with the dust opacity because 
less radiation from the old population is converted into FIR e- 
mission. This smaller influence of the old stellar population on 
the FIR/radio ratio results in a correspondingly smaller effect 
of the long time star formation history on the FIR/radio ratio. 
However, since the UV radiation of massive stars is the domi- 
nant heating source of dust per stellar mass unit (Xu, 1990), we 
expect this effect to be noticeable only for galaxies with very 
low SFR during the past 10 8 yrs. 

The uniform mixing of stars and dust probably leads to a 
decrease of the overall FIR production in a disk galaxy rather 
than to an increase. The reason is that the dominant massive stars 
generally form and evolve in regions of increased dust density. 
This should on average increase the overall galactic optical 
depth for their UV radiation, and thus make galaxies overall 
more "optically-thick" than does uniform mixing (section 5). 

A 1 -dimensional model describes the electron propagation 
in a galaxy correctly if the distribution of the sources of CR's is 
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the disk plays no role (except for regions close to the radial disk 
boundaries) and the electron phase space density, f, is only a 
function of z and E. As long as diffusion is a linear process (no 
effect of the cosmic rays on their own diffusion), this appears a 
reasonable approximation for the integrated radio emission of a 
galaxy. If we were to investigate local correlations between the 
FIR and the radio emissions, clearly a more realistic geometry 
for the distribution of stars and SNe would have to be taken 
into account, and a three-dimensional diffusion model would 
be necessary. 

The synchrotron emission was calculated in a diffusion mo- 
del, not taking into account other modes of CR propagation. 
However, convection can also play a role in the transport of 
CRs, especially at large distances from the disk: For example, 
the observations of the halo radio emission of NGC 463 1 can 
be explained by convection (Hummel & Dettmar 1990). Con- 
vection will be caused by a galactic wind which should exist 
even for a galaxy like ours (Breitschwerdt, McKenzie & Volk 
1991, 1993). It should certainly dominate diffusion in starburst 
galaxies like M82. 

If radiative energy losses dominate the distribution func- 
tion of electrons, the spatial transport mechanism (diffusion 
or convection) affects the integrated spectral index only if the 
magnetic field is spatially inhomogeneous. Otherwise (spatially 
homogeneous B, dominant energy losses), the transport mecha- 
nism is irrelevant for the integrated synchrotron emission. In the 
general case, convection will tend to produce a flatter spectrum 
and a slightly weaker dependence of the synchrotron emission 
on the magnetic field than does purely diffusive transport. This 
can be seen from the following argument: convection is an en- 
ergy independent transport mechanism. Therefore, if the CR 
propagation were only by convection, in the asymptotic case of 
very low frequencies the spectral index a would be the same as 
in Eq.(29) for fi = 

aconv = (^2~^) = °- 6 ( for x = 2 - 2 >- ( 41 ) 

The logarithmic derivative (5 of the synchrotron emission with 
respect to B for \x = is for very low frequencies, 

&onv=^=1.6. (42) 

Thus if convection plays a role, the spectral index will be, also 
at finite frequencies, in even better agreement with the obser- 
vations, because of the flatter asymptotic spectrum of P syn (u). 
The conclusions about the influence of the magnetic field will, 
on the other hand, not be changed significantly. 

6.2. Correlation between [/ ra( j and Ub 

An important result of this paper is the tentative discovery of 
a near proportionality between the energy density of the radi- 
ation field and that of the magnetic field. In turn it could be 
derived from the empirical existence of the FIR/radio correla- 
tion as a necessary condition for it. This additional correlation 
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planation for this correlation could be in terms of turbulent 
dynamo activity in galaxies: The emission from massive stars 
dominates the interstellar radiation field. Through stellar winds 
and SN explosions these stars produce at the same time inter- 
stellar turbulence whose hydrodynamic energy density, Uturb, 
is transported with hydrodynamic velocities that are 10 2 -10 3 
times smaller than the velocity of light. This can compensate 
for the fact that for a quasi-stationary galaxy the time-averaged 
hydrodynamic energy flux into interstellar space is smaller by 
a similar factor than the photon luminosity from the stars. Thus 
the hydrodynamic energy density Uturb ma Y be comparable to 
the energy density U ri d of the interstellar radiation field, as it is 
observed in our Galaxy. In turn, the equipartition of turbulent 
and magnetic pressures limits the dynamo effect that increases 
the magnetic field and therefore can lead to a near equipartition 
between U tur b, U ad and U B - 

A detailed model based on these ideas will be given in a 
future paper. 



7. Conclusions 

We conclude that the FIR/radio correlation is due to a few basic 
properties of late type galaxies like a sufficient optical thick- 
ness of their gas/dust disks, a sufficiently large size of their 
magnetized halos, and a correlation between the energy den- 
sities of the magnetic and the radiation field, most likely due 
to a near equipartition between random kinetic energy density, 
and energy density of the interstellar photon field. Probably the 
disruptive tendency of energy production as a consequence of 
star formation is just in dynamical equilibrium with the ten- 
dency of the gaseous disk and halo to contract gravitationally, 
and to produce stars at a higher rate. Similarly, the energy loss 
of the interstellar medium due to radiative cooling and due to 
the escape of the CR nucleon component (energetically by far 
dominating the radiatively cooling CR electron component) to- 
gether with a galactic wind, just appears large enough to keep 
the various energy densities Uturb , U a &, and U-g, at a comparable 
level. 

In this paper we have only been partly able to explain these 
deep connections. Rather we have theoretically tied together 
their empirical consequences, like the values of optical depths, 
etc.. One of these consequences is the FIR/radio correlation. It 
theoretically requires the existence of an [/ ra d/^B correlation 
which was also found empirically if we accept the minimum 
energy assumption as an approximation to calculate galactic 
magnetic fields. To this extent the U rd d/Us correlation is the 
major discovery of this paper. In principle this correlation al- 
lows the determination of galactic magnetic field strengths from 
observations of the FIR/radio ratio, and of the integrated photon 
flux from a galaxy of known source geometry. In practice many 
other parameters enter this calculation. Nevertheless it is a rare 
case in astrophysics that observations plus theory directly imply 
new astrophysical relations which were not known before but 
are at least vaguely plausible as we have speculated above. 
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nant role of SN precursors and their ultimate explosions for the 
dynamics of the interstellar medium in most normal late-type 
galaxies. However there must be a limit to such a behavior. If 
the current (~ 10 8 yrs) star formation rate becomes too small 
compared to the average rate over times of the order of 10 9 
yrs, then old stars will dominate the FIR radiation without con- 
tributing to the radio continuum emission as argued by Condon 
et al. (1991) and Xu et al. (1994). Thus there is also no phys- 
ical reason any more for a correlation between [/ ra d an d Ub- 
Such galaxies will probably not differ much in their structure 
from the more active ones. But their radiation fields in different 
wavelength ranges will not be correlated, because the energeti- 
cally dominant stellar radiation from old stars does not strongly 
couple to the dynamics of the interstellar medium. From this 
point of view they are less interesting. As a consequence the 
next question regards the opposite case of starburst galaxies 
which are characterized by variations on time scales < 10 8 yrs. 
They also appear to follow the FIR/radio correlation although 
it is not immediately clear why this should be the case. This 
question will be discussed in a companion paper. 
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Appendix A: Solution of the diffusion equation 

In the asymptotic cases \sh\ S> 1 and \sh\ <S 1 Eq. (18) 
is an exact solution to the diffusion equation Eq.(ll). In the 
range between these asymptotes, the function f(E, z) given 
by Eq. (18) does not fulfill Eq. (11) exactly, but is a good 
approximation. Its accuracy was estimated in the following way: 
Eq. (18) satisfies exactly the equation: 

( E \ x 

- 5 z^sn^sn [S(z) + A(£, z)] , 

\m e c z J 

where the residue A(E, z) equals: 

A(E,z) = 2^ 1 F 1 (p,\, s) 

( \F\(y+ j, |,s fe ) giEKp, j,s h ) 
V h i^ifei^) 2 ds h 

1 iF l (p+ \, \,s h ) _ s h diFi(p+ j,j,s h ) \ 

2 ,Fi(p, \,s h ) x F x (p,\,s h ) ds h J 




log(i/)[MHz] 



Fig. Al. The integral over the residue J Zh A(E, z)dz. The parameters 
describing the halo geometry are the same as in Fig. 3, i.e. L/zh =0.1, 
5 = p = 1 (full line), 5 = 0.5, p = 0.4 (dashed line) and 5 = 0.2, p = 0.4 
(dashed-dotted line). 

For the case of S < 1 , or p < 1 , an analogous (but more compli- 
cated) expression for A(E, z) can be derived. The expression: 



f_ h Zh HE{v) jZ )Az 
C Zh S(z)dz 



r 

J -z h 



A(E(p),z)dz 



gives an estimate for the error of the approximate solution. s(v) 
is shown in Fig. (Al), both for the spatially homogeneous as 
well as for the inhomogeneous case. It can be seen that the 
accuracy in the homogeneous case is always better than 12 %, 
and in the inhomogeneous case it is better than 14 %. The error 
of a syn and of the derivatives of Qia9GHz listed in Table 1 is 
much lower: it is for all quantities less than 2 %. 

Appendix B: Dust opacity 

With the data available for the galaxy sample, we want to cal- 
culate the dust opacity and its dispersion. The sample of 1 14 
normal spiral galaxies that we use in this paper is the same as 
used and described in Xu et al. (1994). It is selected from the 
catalogue of Aaronson et al. (1982) of nearby late-type galaxies 
with known H-band (1.6/um) magnitude, for which also FIR 
and radio fluxes are available. 

We are interested in the diffuse dust, associated with HI, 
which in our Galaxy has been shown to emit a large fraction 
of the total FIR emission (Boulanger & Perault 1988, Cox & 
Mezger 1988, Perault et al. 1988, Bloemen, Deul & Thaddeus 
1990). The FIR emission itself is a poor indicator of the dust 
mass because it depends very sensitively on the, not very well 
known, dust temperature distribution. More importantly, it does 
not allow a distinction regarding the origin of the dust emission 
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approach: We assume that the dust-to-gas ratio, which for our 
Galaxy is observed to be quite constant (e.g. Savage & Mathis 
1979), is the same in other galaxies. For some nearby galaxies 
this assumption has been shown to be valid: Rice et al. (1990) 
could fit the infrared emission and the optical and UV spectrum 
of M33, using the Galactic value of the dust-to-gas ratio. For 
M3 1 , Walterbos & Kennicutt (1988) showed that in M3 1 the val- 
ue is also the same as in our Galaxy. In the very metal poor Large 
Magellanic Cloud, the ratio is higher by a factor of 4 (Koornneef 
1982) however, and in the Small Magellanic Cloud even higher 
by a factor of 8-12 (Bouchet et al., 1985). Since the galaxy 
sample used here consists mostly of spiral galaxies, we believe 
that a constant gas-to-dust ratio is a reasonable assumption in 
this case. We take N(HI)/E(B - V) = 5 • 10 21 atoms/cm 2 mag 
(Savage & Mathis 1979) with N(HI) being the column densi- 
ty of HI, and E(B-V) the color excess. This gives an average 
absorption cross section per H-atom: 



R E(B - V) 
1.086 N(HI) 



= 5.71 • 10- 22 cm7atom 



(Bl) 



with R = 3.1 (Savage & Mathis 1979). The visual dust opacity 
rv in a galaxy can therefore be estimated with the use of cry, 
the HI mass F(HI), and the diameter of the galaxy D: 



Ty = ay- 



F(HI) 



1 



m p tt(D/2) 2 



(B2) 



m p is the proton mass. The opacity in the UV is higher: The 
average extinction curve of dust (Savage & Mathis 1979) shows 
that the extinction at 912 A is a factor of 5.72 higher than 
in the visual wavelength range. For the galactic diameter we 
can use the optical isophotal diameter D25 instead of the HI 
diameter with a small correction factor: Fouque (1983) found 
that the effective HI diameter D H , i.e. the diameter within which 
half of the HI mass is situated, and D25 are nearly the same: 
The ratio D^/Dis = 1 — 1.5 for spiral galaxies, depending on 
morphological type, with a very low dispersion of less than 0. 1 
in the logarithm. Thus we approximate the FIR disk diameter 
by D25 and apply the morphological type dependent correction 
factor of Fouque. 

The average opacity determined in this way at 2000 A 
is < tuv >= 1-1 with a dispersion in the logarithm of 
<5iogr = 0.27. There is no tendency noted for the opacity to 
increase or decrease with luminosity. Thus the galaxies are 
optically thick for the non-ionizing UV light which is mainly 
emitted by massive stars. In the optical wavelength range galax- 
ies are optically thin (< ry >= 0.4). Therefore preferentially 
the radiation from massive stars is absorbed. 

Appendix C: Calculation of the energy density of the radi- 
ation field and of the magnetic field for the galaxy sample 

Assuming that the radiation intensity in the disk is homogeneous 
and isotropic, the energy density of the radiation field in the 
galactic disk is given by: 



Utad = 



J bol 



C7r(£> 25 /2) 2 (D 25 /kpc)' 



•5.56-10 



-35 



(CI) 



with D25 being the optical isophotal diameter at 25 mag, and 
L bo i the bolometric luminosity. The photons contributing to the 
inverse Compton losses have to be below the Klein-Nishina 
limit which means that they have to follow the condition: 



hv E 



< 1 



(e.g. Longair 1992). For electrons in the energy regime con- 
sidered here (about 1—20 GeV) the upper limit lies in the 
X-ray range (hv ss 200 eV). Data for the luminosities at vari- 
ous wavelength bands are available: the blue and H magnitudes 
and the FIR emission. With these data, and taking into account 
the microwave background, the total radiation energy density 
of a galaxy can be estimated by extrapolation. From the blue 
magnitude the magnitudes in the visual and in the UV are esti- 
mated by the average (U-B) and (B-V) values corresponding to 
each morphological type (Huchra 1977). Obtaining in this way 
the luminosities per wavelength at UV, visual, blue wavelength 
and at 1.6 /xm (H-band), one can estimate the total luminosi- 
ty between 912 A and 4 /tm by extrapolation, referring to the 
calculated galaxy spectra in Mazzei, Xu & De Zotti (1992). 



£912-4000 [W] = L uv [W/nm] ■ Q.lfim + L B [W/fj,m]- 
0.1/xm + L v [W/iim] ■ 0.1/xm + L H [W/ /im] • 4/tm. 



(C2) 



In addition we take into account the FIR luminosity (40 — 
120/im) and the microwave background with an energy density 
of f/back = 0.25eV/cm 3 . The wavelength range between 4/tm 
and 40/um is negligible for the total luminosity of a galaxy. 
The total radiation energy density is the sum of the energy 
densities at different wavelengths. All these energy densities 
individually correlate positively that the energy density of the 
magnetic field, but with very different slopes (between 0.6 and 
1.7). As a consequence, the slope in the correlation between 
[Trad and Ub depends rather sensitively on the weight which 
each luminosity has in Eq. (Bl). The correlation coefficient 
and the dispersion of the correlation, on the other hand, are not 
significantly affected by it. 

The magnetic field is estimated by the standard minimum 
energy approximation (e.g. Pacholczyk 1970). It is based on 
the assumption that the value of the magnetic field is such that 
the sum of the CR energy and the magnetic field energy is 
minimal. The magnetic field can in this way be calculated from 
the synchrotron emission P syn (v) and the volume of the radio 
emitting region: 



= 4.2 10 



(P 5yn (i/)/W/Hz)(i//GHz) a 
(D/kpc) 2 



[HGI 



In this formula we used the standard parameters: The propor- 
tionality factor between the energy of the cosmic ray electrons 
and the total energy of the cosmic rays is k = 1 00, the spectral in- 
dex of the synchrotron emission a = 0.7 and the minimum and 
maximum frequency of the synchrotron emission v min = 10 
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to be a cylinder of 1 kpc height. Its diameter D is given by the 
lowest contour line of the galaxy maps given in Condon (1987) 
and Condon, Yin & Burstein (1987). Only 86 galaxies of the 
sample could be considered because for the rest no diameter 
could be deduced. 
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